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We investigate the possibility of the electrical control of spin transfer in monoatomic chains 
incorporating spin-impurities. Our theoretical framework is the mixed quantum-classical (Ehrenfest) 
description of the spin dynamics, in the spirit of the s-d-model, where the itinerant electrons are 
described by a tight-binding model while localized spins are treated classically. Our main focus is 
on the dynamical exchange interaction between two well-separated spins. This can be quantified by 
the transfer of excitations in the form of transverse spin oscillations. We systematically study the 
effect of an electrostatic gate bias V g on the interconnecting channel and we map out the long-range 
dynamical spin transfer as a function of V g . We identify regions of V g giving rise to significant 
amplification of the spin transmission at low frequencies and relate this to the electronic structure 
of the channel. 
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I. INTRODUCTION 

The rapid development of the field of spintronics^ over 
the past two decades has uncovered exciting and novel 
phenomena related to the dynamics of the electronic spin 
in a wide variety of systems, ranging from bulk materi- 
als to spatially-confined structures^. Fueled by the ever- 
growing needs for speed, capacity and energy-efficiency in 
computing, the central objective in understanding and ul- 
timately controlling the spin properties in the solid state 
has been constantly shifting towards the nano-scale. At 
these lengths and times the conventional methods for spin 
control, based on magnetic fields alone, are limited by 
scalability issues. Alternative approaches are thus sought 
and they typically involve electric fields of some forrrP. 

One way to manipulate spins by purely electrical 
means relies on the spin-transfer torque mechanism 4 . 
This approach uses spin-polarized currents to control the 
direction of the magnetization and has been realized in 
various nano-structured materials ranging from magnetic 
multilayers 5 to, more recently, single atoms in STM-type 
geometries^. Alternative to electric current control is the 
optical control, such as in the laser-driven ultrafast mag- 
netization switching^. In a somewhat different context, 
the optical manipulation of single spins in bulk media 9 
is at the heart of the most promising candidates for the 
quantum information processing technology * 10 * 11 !. 

Another alternative is based on the idea of an electro- 
static control of the spin-density, i.e. of the construction 
of spin-transistor type devices^. Recently, the concept 
of gate-modulated spin-pumping^ transistors has been 
studied theoretically in infinite graphene strips with pat- 
terned magnetic implantations 14 . Importantly, such de- 
vices rely on the efficient transport of spin information 
between two points in space and time, and require the 
possibility to actively tune the propagating spin-signal 
during its transport. In this work we explore this possi- 
bility for atomistic spin-conductors. We consider a finite 
mono-atomic wire linking two localized spin-carrying im- 
purities. When one of the localized spins is set into pre- 



cession, it generates a perturbation in the spin-density. 
This perturbation is carried through the wire by con- 
ducting electrons and can be detected in the dynamical 
response of the second spin. We show that the propaga- 
tion of the spin-signal and consequentially the dynami- 
cal communication between the two spin centers, can be 
tuned by means of an electrostatic gate applied to the in- 
terconnecting wire. Our main finding is an enhancement 
of the communication for a certain range of gate volt- 
ages. This is linked to the modification of the electronic 
structure of the wire induced by the applied electrostatic 
gate. 

Because of their considerable size the systems investi- 
gated here are still beyond the present numerical capabil- 
ities of first-principles dynamics^ and as such they are 
described by model Hamiltonians. Usually the dynam- 
ics is approached within the linear response approxima- 
tion^. In this paper we go beyond the linear response 
limit and propose a fully microscopic description based 
on the time-dependent Schrodinger equation, which al- 
lows us to describe arbitrary excitations. In other words 
our simulations are not limited to small gate voltages or 
small-angle spin precession. We use a single-band tight- 
binding Hamiltonian to model the itinerant s-electrons 
in our metallic wires and include local Heisenberg inter- 
actions to a number of magnetic ions (typically one or 
two). The latter are modeled as classical spin degrees of 
freedom and enter on equal footing in the common mixed 
quantum-classical dynamic portrait of the system. 

This paper is organized as follows. In Section [ll] we 
introduce the model system and our theoretical frame- 
work. Section |III| contains the results of our investiga- 
tions. Firstly we investigate the electron response of an 
atomic wire, not including magnetic impurities, to lo- 
cal spin excitations. Secondly, we study the dynamical 
interplay between two spin-impurities electronically con- 
nected by the wire and show how such a dynamics is 
affected by the gate voltage. Finally we draw some con- 
clusions. 



II. THE MODEL 

A cartoon of the device considered is presented in Fig. 
[I] This consists of an TV-sites long atomic wire intercon- 
necting two magnetic centers. The latter, represented in 
the figure by the spin- vectors, enter our model as substi- 
tutional magnetic impurities positioned at the two ends 
of the wire (the local spin originates from the deeply lo- 
calized d-electrons). One of them, say the left-hand side 
spin-center Si , is labeled as "driven" , as its precession is 
induced and sustained by a local (at that site) magnetic 
field. The other localized spin, £2, located at the other 
end of the chain is the "probe" spin and it is not directly 
coupled to any external magnetic field but only to the 
electron gas. In this way £2 probes the magnetic excita- 
tions produced by the driven spin Si as these propagate 
through the interconnecting wire. The device is described 
as a mixed quantum-classical system in the spirit of the 
(s-d)-modeP^, where the time-dependent Hamiltonians 
of the two exchange-coupled spin sub-systems read 
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FIG. 1: (Color online) Model system investigated in this work: 
two localized spins are electronically connected by a mono- 
atomic wire, where electrons can flow. An electrostatic gate 
is applied to some of the interconnecting sites. 
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The top expression is for the quantum electrons. Here 
Hj 3 B = SiSij +7<Si,i±i is a single-orbital tight-binding 
Hamiltonian with on-site energies S{ and hopping inte- 
gral 7 (7 sets the relevant energy scale for the entire sys- 
tem); (cf) is the creation (annihilation) operator for 
an electron with spin-up (a = 1) or spin-down (a = 2) at 
the atomic site i; a = ^(cr x ^a y ^a z ) is the electron spin 
operator, {cri}i= x , y ,z being the set of Pauli matrices; Si^ 
is the unit vector in the direction of the localized spin; 
J > is the exchange coupling strength. 

The classical dynamics of the local spins is governed 
by Hs(t), which describes the interaction of Si^ with the 
mean-field local electron spin-density Si = (<r) i5 taken as 
the instantaneous expectation value of the conduction- 
electron spin at site i (see further in text for the exact 
definition). The classical Hamiltonian also includes in- 
teraction with the external magnetic field B = (0, 0, B z ), 
which is applied locally to Si only, in order to drive its 
precession. We further assume g = 2 for the localized 
spins and [1b = 5.788 10 _5 eV/T is the Bohr magneton. 

In order to study the dynamics of this system in the 
time-domain we integrate the coupled quantum and clas- 
sical e quatio ns of motion (EOM) of the two spin sub- 
systems^CL^. In our spin-analogue to Ehrenfest molecu- 
lar dynamics^, the full set of coupled Liouville equations 



Here {, } represents the classical Poisson bracket, [, ] 
stands for the quantum-mechanical commutator and S = 
h is the magnitude of the two classical spins. The first 
EOM is for the electron density matrix p. At the initial 
time, to , p is constructed from the eigenstates {|^)}^=i 
of the spin-polarized electron Hamiltonian H e \(to) [see 
Eq. ([2])], with the composite index v = {i,cr} label- 
ing the set of 2N spin-polarized eigenstates. We define 
P(to) = fA { Pv){Vv\ with f v = r] F (e u - E F ) being 
the occupation numbers distributed according to Fermi- 
Dirac statistics. The instantaneous onsite spin-density is 
thus generated as 

Si (t) ee (t) = Tr [p(t)a] t = £ pf{t)a^ . (4) 

The coupled EOMs are integrated numerically by using 
the fourth-order Runge-Kutta (RK4) algorithm 20 . As a 
result the set of trajectories for the local electronic spin- 
densities Si(t) are obtained as well as those of the driven 
and the probe classical spins Si^(t). The typical length 
of the chain that we consider is N = 100. 

The effect of an electrostatic gate is incorporated in our 
model as a rigid shift of the onsite energies Si —> Si + V g , 
where i G [^1,^2] is a certain range of sites in the middle 
of the chain and V g is the gate voltage^. 
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III. DYNAMICS OF THE ITINERANT SPINS 
WITH FROZEN IMPURITIES 



sub- femtosecond duration) at least over the investigated 
time-scales of a few picoseconds. 



A. No external gate 

As a preliminary step towards the combined quantum- 
classical dynamics we first address the dynamics of the 
spin-density of the itinerant electrons in the presence 
of the two local spins, whose directions are fixed, e.g. 
Si^H^- A spin excitation is produced by a small but fi- 
nite spatially-localized perturbation in the spin-density. 
As initial density matrix at to we use the one that corre- 
sponds to a perturbed Hamiltonian H e \(to — St) in which 
one of the localized spins is slightly tilted in the x — z 
plane, such that Si (to — St) • z ~ dO < 5°. However at 
to we bring Si back to its original direction (\\z) where 
it stays throughout the simulation. In other words we 
study the time evolution of the system with the two lo- 
cal spins parallel to each other starting from the ground 
state electronic charge density of the system where one 
of the two spins is tilted by a small angle. 

The evolution of the density matrix for t > to is then 
given by p(t) = e~ tHelt / h p(to)e %Helt / h , which translates 
into the following expression for the density matrix ele- 
ments 
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(5) 

In Eq. (I5| we have defined the frequencies uj mn = e m — 
e n )/h corresponding to differences between the eigen- 
values e m of H e \(to). The coefficients c£ m = (ka\ tp m ) 
are the projections of the eigenvectors \(p m ) on the spin- 
resolved atomic orbital basis \kcr) = \k) <g) |cr), where k 
represents the atomic site and a =t, 4- is the spin compo- 
nent. 

The typical evolution of the spin density of the itiner- 
ant electrons resulting from excitation described above is 
presented in Fig. [2] The trajectories of the individual on- 
site spin polarizations Asf (t) = sf (t) — sf (to) stemming 
from from Eq. ^ are perfectly identical (on the scale 
of the graph) to those obtained by the numerical inte- 
gration of the EOM [i.e. from Eq. ([2|], confirming the 
reliability of our time-integrator for the typical duration 
of the simulations. 

The spatial distribution of the initial spin polarization, 



sf(to), is shown in Fig. 
tion, presented in Fig. \2. 



2Tb) and its subsequent evolu- 
^c) can be qualitatively charac- 
terized as the propagation of a spatially-localized spin 
wave-packet. This travels along the wire with a practi- 
cally uniform velocity very close to the Fermi level group 
velocity, as expected in view of the rather small over- 
all local spin-polarization of the chain. The packet then 
gradually loses its sharpness as it disperses. However, 
the backbone of the packet is detectable even after a few 
tens of reflections at the ends of the wire. Such a feature 
demonstrates that this finite ID model atomic system is 
a relatively efficient waveguide for spin wave-packets (of 
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FIG. 2: (Color online) Spin excitation of one dimensional 
wire, (a) Time dependence of the spin-polarization on the 
first site (Asf) as obtained from the numerical integration 
of the EOMs and directly from Eq. |5]). (b) Initial spin- 
polarization s z (to) as a function of site index, (c) Time and 
space evolution of Asf with the color shade representing the 
magnitude of absolute value of Asf(t). Note that 2r is the 
wave-packet round-trip time as also seen in the bottom panel. 

Eq. (§ can also be used to calculate directly the spec- 
tra of particular spin observables. These match perfectly 
with those calculated by performing the discrete Fourier 
Transform^ (dFT) of the time-dependent spin-densities 
obtained over the finite duration of the dynamic simu- 
lation. Clearly, in the absence of the localized impuri- 
ties (i.e. for a finite homogeneous tight-binding chain) 
uj mn can be calculated exactly. The spectrum of Asf 
has monotonously decreasing amplitudes from the lowest 

possible frequency u; m i n « p^pjy^ to the maximum one 

w max ~ -^p (expressed in the limit of TV — >> oo). This is 
also true for the case of interaction with the two frozen 
local impurities (e.g. for J = 7) as these have a minor 
effect on the electron Hamiltonian. For the parameters 
typically used in our simulations the corresponding max- 
imum period T max = -^ JL — < 1 ps is within the total time 
of the simulation while the corresponding minimum pe- 
riod T m i n = ^ 27r w 1 fs is much larger than the typical 
time step At = 0.01 fs. 

A more detailed spectral analysis is obtained by the 
two-dimensional dFT power portraits 21 of As£ (see Fig. 
|3| denoted as dFT[Asf (t)](/c, u). In Fig. |] we com- 
pare such spectra for the itinerant spin-dynamics to its 
exact counterpart in the case of a finite homogeneous 
tight-binding chain. These portraits reveal key features 
of the one-dimensional fermionic system that vary sys- 
tematically with the band-filling po? namely (i) a near- 
continuum of allowed modes in a certain (fc, cj)-space 
region, defined by a low- and a high-energy dispersion 
functions, (ii) a linear dispersion for small k (uo ex k for 
k —> 0). An analogous mode-occupation patterns have 
been rigorously analyzed in relation to the dynamical 
properties of one-dimensional quantum Heisenberg spin 
chains which too have a tight-binding type dispersion re- 
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latiorPH The low-energy mode-occupation limit is due 
to the fact that the energy of the electron-hole excita- 
tion can approach zero only for Ak — >• and Ak — >• 2fcp, 
where /cp is the Fermi vector (fcp = ^ f° r Po = 0.5, cor- 
responding to one electron per site). The variation of the 
band filling away from the half-filling results into a fold- 
ing of the low-energy limit as shown in Figs. [3^b), [3^c), 
|3|e) and|3|f). Note that due to electron- hole symmetry 
we only show the spectra for po > 0.5. 
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FIG. 3: (Color online) dFT[Asf (t)} (k, u) for different band- 
filling: (a, d) po = 0.5 , (b, e) po = 0.6 and (c, f) po = 0.8. 
The top panels show the exact analytical excitation spectrum 
for a homogeneous tight-binding chain without local spins. 
The bottom panels show the results of our numerical sim- 
ulations. The inter-site distance a is arbitrary. The color 
shade in the bottom panel represents the absolute value of 
dFT[Asf (£)](&, oS) and is logarithmically scaled for better con- 
trast. 



B. Electrostatic gate applied 

The dynamics of the itinerant spins, as described by 
Eq. ([5]), in the case of an electrostatic gate applied to a 
section of the wire (typically in the middle of the wire) 
cannot be expressed in a closed form as a function of 
V g for arbitrarily big systems. As such we resort to our 
numerical integration scheme. Before addressing the dy- 
namics, however, we first analyze the ground-state elec- 
tronic structure of the gated wires for different values of 
the gate potential V g . Displayed in Fig.[4]is the adiabatic 
variation of the eigenvalues for a non spin-polarized chain 
with N = 30 sites (without localized spins) as a function 
of the gate potential V g (the gate is applied to 10 sites 
in the middle of the chain, i.e. at the sites with indexes 
going from ii = 11 to 12 = 20). Note that energy-related 
quantities on all the figures are in units of the hopping 
integral 7. For V g close to 0, the discrete energy spec- 
trum spans in the range [—27,27]. With the increase of 
V g the levels spacing distorts as the eigenstates are af- 
fected differently by the gate. Certain eigenvalues grow 



nearly linearly with increasing V g (these are shown as red 
squares in Fig.[4|. The spatial distribution of these eigen- 
states is predominantly concentrated in the gated region, 
i.e. they corresponds to the region where the local onsite 
energy has been modified. 
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FIG. 4: (Color online) Schematic of the gate-dependence 
of the ground state energy spectrum of H e \. According to 
their weights at the three subsections of the chain Qc,m,k — 
^2 ieC M n \cin\ 2 we distinguish three types of states depend- 
ing on which of the three partial weights is the greatest. We 
use black circles for the case of Qc, red squares for Qm and 
green diamonds for f^. Note that the eigenstates correspond- 
ing to the two ungated regions (Qc and Q-r,) are degenerate 
by symmetry, as the gate is applied in the exact center of the 
wire. The inset shows a magnification of an area of intermedi- 
ate V g , illustrating the situation of avoided crossings. A short 
chain with N — 30 sites is used for simplicity. 

For extremely large values of V g (V g > 47) the chain is 
effectively split into three energetically-decoupled parts, 
the gated middle and the two identical un-biased ends on 
the left-hand side and on the right-hand side. More inter- 
esting for us is the range of intermediate gate voltages, 
for which the three parts of the wire are substantially 
affected but not yet decoupled by the gate voltage. For 
such V g avoided crossings occur between states localized 
in the gated and non-gated regions. This gives rise to 
additional low-frequency lines in the dynamical spectrum 
as described by Eq. ([5|. As we will demonstrate in the 
Appendix the presence of these avoided crossings around 
the Fermi level for certain intermediate V g yields an en- 
hanced transmission through the gated wire (waveguide) 
at low frequencies. 

The two-dimensional dFT images (Fig. [5| of the spin 
dynamics in the presence of the gate bring additional di- 
mension to the electron spectroscopy analysis. The dif- 
ference here with respect to the case depicted in Fig.|3|d) 
is that a gate has been applied to the middle of the chain 
in the ground state, i.e. V g i2 ] c% has been 

added to H e \ at t = to- Again we use the half-filling case, 
po = 0.5 (one electron per atom). From the figure it is 
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immediately noticeable that the effect of the variation of 
V g on the excitation spectra is quite similar to the effect 
of the band- filling in the non-gated case. Due to the gate 
potential, the relative electron populations of the gated 
and gate-free parts of the chain change (the gated re- 
gion is depopulated). For intermediate values of V g [see 
Fig. |5jb)] the excitation spectrum is rather a superpo- 
sition of two spectra with band fillings below and above 
0.5 (approximately 0.2 and 0.6). Furthermore we find a 
substantially increased population of the low-frequency 
modes for all /c-vectors, i.e. of states were forbidden by 
symmetry for V g =0. For large V g [see Fig. 5jc)] the 
middle part of the wire becomes almost completely de- 
pleted and the (k, lj) -port rait corresponds effectively to 
the excitation spectrum of a single chain with a band- 
filling po w 0.75, which is similar to the non-gated case 
presented in Fig. [3F 
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dFT images for different values of the gate voltage are 
presented in Fig. [6j 
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FIG. 5: (Color online) dFT[Asf (t)](k, u) for different values 
of the gate potential: (a) V g =0.27 , (b) V g = 2.2^ and (c) 
V g = 4.67. The color shade is identical to that on Fig. [3] 



IV. COMBINED QUANTUM-CLASSICAL SPIN 
DYNAMICS IN THE PRESENCE OF A GATE 

The inclusion of dynamic local spin-impurities in the 
model requires the numerical integration of the set of cou- 
pled non-linear EOMs of Eq. ([3|. For a small number of 
classical spins the dynamics of the itinerant spin-density 
is qualitatively very similar to what is described by Eq. 
([5]). The frequency- domain analysis of the classical spin 
trajectories by means of dFT reveals the characteristic 
signature of the discrete electronic spectrum with only 
additional modulations in the amplitudes. We focus now 
on the case in which Si is driven by a local magnetic field 
into a precession, hence it acts as a spin-pump (see car- 
toon in Fig.[T]). It is important to note that what we refer 
to as a local magnetic field B = (0, 0, B z ) is only instru- 
mental to trigger and sustain a uniform Larmor preces- 
sion of Si, i.e. it does not produce any Zeeman splitting 
in the itinerant electrons spectrum. We consider now the 
following situation: the quantum-classical spin-system is 
in its ground state until t = to, when the first local spin 
Si starts fluctuating to form a small misalignment with 
B. This sets the entire quantum-classical spin system 
into motion. The typical trajectories of the transverse 
components Sf(t) and Sf (t) and their frequency-domain 
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FIG. 6: (Color online) Time evolution of the x-components of 
the localized spins, Sf (t) (blue) and Sf (t) (red) [left] and the 
corresponding spectra, dFT[Sf (t)](cj), [right] for three differ- 
ent gate voltages (a, d) V g = O.I7, (b, e) V g = 2.2j and (c, 
f) V g = 4.67. The frequency- domain is represented by the di- 
mensionless quantity fiw/'y and ujl = huuL/'J = 1.16 for 7 = 1 
eV. Note that in the case of V g — 2.2j different scales are 
used for the magnitudes of Sf(t) and Sf (t) and their spectra 
(axes corresponding to Sf (£) are marked in red). 



In this case the excitation of a transverse spin-density 
at t — by the initiation of the Larmor precession is very 
similar in nature to the excitation ind uced by tilting one 
local spin investigated before (section |TlT| ). In fact it in- 
deed develops in a very similar way in the early stages of 
the time evolution. Just like in Fig. [2] a non-equilibrium 
transverse spin-density spin packet is sent along the wire 
and both the localized spins respond each time the packet 
reaches them. Figure |6ja) displays the case of a very 
small gate voltage V g <C 7. The frequency-domain image 
of the precessional motion of Sf(t) has the signature of 
the discrete electronic spectrum with an amplitude enve- 
lope that peaks at the Larmor frequency cjl = 2/i#£?i/fL 
The dFT-portrait of the probe spin Sf (t) is very simi- 
lar to that of Sf(t). Although significantly down-scaled 
in amplitude, it has the electronic-structure modes im- 
printed in its classical motion in the same way as it is 
for the driven spin. This qualitative behavior is observed 
for V g = and a range of small voltages V g < 7 but 
it changes dramatically for gate voltages V g > 7 [see 
Fig. [6jb) and (c)]. The case of extremely big voltages 
V g > 47 (depicted in the bottom panel) is a trivial one for 
which the electrostatic barrier V g is completely opaque to 
the electrons and the probe spin S2 does not respond to 
the Larmor precession of Si (note also the rarefaction of 
modes in the Sf spectrum due to the effective shortening 
of the chain by 2/3). 

Special attention must be devoted to the intermediate 



Vg regime in which the dynamics of the probe spin is qual- 
itatively different and uncorrected to that of the driven 
one [see Fig. |6jb)]. The typical outcome of the real time 
dynamics at such intermediate voltages 7 < V g < 47 is 
that the probe spin starts accumulating very big trans- 
verse deflections. This is then fed back into the pumping 
spin which also deflects more but still preserves, to a 
great extent, its precession about the magnetic field. 

It is worth noting that this self-amplification of the 
spin-dynamics does not come at a total energy cost, as 
the model system described by Eg. Q is completely 
conservative. As such we simply observe a conversion 
of electrostatic energy into "spin-energy" in the form 
of a spin amplitude transverse to the driving magnetic 
field. This accumulation is related to the increasing en- 
ergy transfer from the itinerant electrons to the localized 
moments, as illustrated in Fig. (7^a). Such an energy 
transfer is defined as AE(t) = \AE e \(t) — AEs(t)\, where 
AE el/s (t) = E el/S (t)-E el/S (t ) : E el (t) = Ti[p(t) • H el (t)] 
and Es(t) = —gHBS\{t) • B. The total energy conser- 
vation AE e i(t) = —AEs(t) has been verified within a 
relative error of 10 -9 %. The amplification of the energy 
transfer^ is characteristic of this V g range (7 < V g < 47) 
[see the inset of Fig. (7^a)]. 

The frequency-domain image of the localized spin tra- 
jectories [see Fig. [6Fe)] for 7 < V g < 47 shows a substan- 
tial qualitative difference for the driven and the probe 
spins. While around cjl the shape of the spectrum of 
Si is similar to the two extreme V g cases, a significant 
amplitude is accumulated now at low frequencies. The 
main difference from the small V g case, however, is in the 
spectrum of £2 which shows practically no response at 
the Larmor frequency but has very large amplitudes in 
the low frequency range. 

This behaviour in the classical spin dynamics charac- 
teristic of certain gate voltages stems from the dominance 
of the low frequency modes in the electron dynamics at 
those V g . The effect can be seen also in the much simpler 
case of a non-spin-polarized uniform atomic wire sub- 
jected to a sudden switch of an electrostatic gate in the 
middle. In the Appendix we consider the charge dynam- 
ics associated to this situation and analyze the lowest 
frequency components entering in Eq. The initial 

charge excitation, similarly to the transverse spin polar- 
ization in the model above, is triggered by a local elec- 
trostatic perturbation at the first site. It follows from 
Eq. ([5]) that, for certain values of V g for which two adja- 
cent adiabatic eigenenergies (e^e^+i) around the Fermi 
level come closer to each other in an avoided crossing, the 
dynamical charge amplitude at the opposite end of the 
wire exhibits a peak at that particular very low frequency 
(e i+ i — 6i)/H. In other words, a certain gate voltage con- 
dition is created for resonant charge transfer between the 
three parts of the chain at low frequencies. This, in the 
case of itinerant transverse spin-polarization, promotes 
the enhanced reaction of the probe spin to the spin pump- 
ing produced on the other side of the gate. 

The electronic part of full quantum-classical spin dy- 



0.09 
^0.06 
< 0.03 



(a) $ - 1 



T 1 1 1 1 1 1 1 1 — 



V=2.3y 



1 2 3 4 5 

I . 1_ ^AlrfLli* 

500 1000 
t(fs) 




1500 



2000 



Q 0.6 
!s 0.4 
°- 0.2 



(b) 



— left 

— middle 
ght 



500 



1000 

t(fs) 



1500 



2000 



FIG. 7: (Color online) (a) Time evolution of the energy 
transfer, AE(t), between electronic and localized spins for 
V g — 2.27. The inset shows the maximum value of the energy 
transfer, AE max , for a fixed duration of simulation (2ps) as 
a function of V g . A red arrow indicates the particular voltage 
for which AE(t) is plotted, (b) Evolution of the electron pop- 
ulations of different parts of the chain for V g = 2.2 7. The two 
black dashed lines show the equilibrium populations of the 
gate-free parts (upper line) and of the gated middle section 
(lower line), i.e. the populations corresponding to the case 
when the gate is applied in the ground state. 



namics in the presence of the gate (Fig. |6| is, in many 
respects, similar to the charge-only case, described in Sec- 



tion IIIB and in the Appendix. As before we construct 
the two-dimensional dFT portraits of the spin-density 
evolution. Providing additional /c-resolved information 
these give another perspective for interpreting the dra- 
matic increase of the low-frequency oscillations of the 
probe spin at certain gate voltages, when compared to 
the case of frozen Si and electrons at t = to — 5t relaxed 
to the external gate potential (Fig. [5]). In the latter case, 
the increased low-frequency occupations (with respect to 
the non-gated case) were attributed to the split of the 
one-dimensional system into subsystems with different 
average electron densities. Here we find a similar dFT 
pattern, now in the presence of dynamically-coupled lo- 
cal spins and a gate abruptly introduced at t = to . Those 
two additional time-dependent factors result in a signif- 
icant variation in time of the electron distributions in 
the gated and non-gated segments of the chain [see Fig. 
[Tjb)]. The dFT image in this case can be interpreted 
as a superposition of spectra corresponding to different 
po, resulting in a smearing of the low- frequency structure 
previously present in the (fc, cj)-space image. A compar- 
ison of the results for V g = and V g = 2.27 is presented 
in Fig. [8] We find that the full dynamical V g = por- 
trait is rather similar to the results obtained for frozen- 
impurities (see Fig. [5| , the only additional feature being 
the adiabatic excitation at the Larmor frequency of the 
driven spin. The difference (point-by-point subtraction of 
the contour plots) with the corresponding gated case, de- 
picted Fig. [8^b) shows that it is indeed mainly the lower 
frequency band of the spectrum (for any wave-vector) 
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that gains occupation in the gated case. 




ka/7i 



FIG. 8: (a) dFT[Asf (*)](&, u) for V g = 0. (b) The 
difference between dFT[Asf cu) for V g = 2.2j and 

for V g = defined as A = |dFT[s?(t)](fc,cj)v fl = 2 .2 7 - 
dFT[sf (t)](k, Lj)v g =o\- The color shade represents the abso- 
lute value of dFT[Asf (£)](&, cj) and A, respectively, and is 
logarithmically scaled for better contrast. 



V. CONCLUSIONS 

In summary, we have employed atomistic dynamical 
simulations to investigate the effect of an electrostatic 
gate as means of controlling the indirect coupling be- 
tween two distant localized spin impurities in a finite 
metallic wire comprising one hundred atoms. One of the 
impurities, precessing in external magnetic field, plays 
the role of a spin-pump and the response of the second 
(probe) spin was analyzed as a function of the gate po- 
tential applied to the interconnecting wire. We identified 
a range of external gate potentials for which the spin- 
pumping is extremely efficient and leads to a substantial 
excitation of transversely-polarized itinerant spin density 
in the non-gated parts of the chain. This, in turn, pro- 
duces a massive low-frequency swing of the probe spin. 
Such a resonant effect has been related to gate-induced 
avoided crossings in the electronic structure of the in- 
terconnect. For certain gate voltages these occur in the 
vicinity of the Fermi level and assist the enhanced trans- 
fer of charge/spin across the barrier. Evidences for this 
effect are also identified in the Fourier portraits of the 
calculated time-dependent spin-distribution. 
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Appendix 

We consider the eigenstate-resolved dynamical ampli- 
tude A mn corresponding to a mode with frequency cj mn , 
associated to the transitions between the eigenstates e m 
and e n . The system is the simplified wire used for the 
spectroscopic analysis of Fig. [4j i.e. it is a non-spin- 
polarized 30-atom long chain. Analogously to the spin- 
polarized case, the t = to excitation is applied simply as 
a potential shift at the first atomic site. From Eq. ([5| 
the dynamically-driven charge- accumulation amplitude 
at site i is defined as 

A mn = c im c in^2 C tm C lnPkl(to) , (A.l) 
kl 

such that 

fti (*) = £4 l ne- <&W (A.2) 

ran 

is the dynamic electron occupation at site i. This quan- 
tity is calculated at the last site, i = iV, i.e. at the oppo- 
site side o f the chain with respect to the perturbation. In 
Eq. ( |A.1[ ) Cim = (i\(pm) are the onsite components of the 
eigenvectors \<Pm) of the non-spin-polarized Hamiltonian 
H e \(V g ) for which the gate voltage V g is applied (between 
the sites i = 11 and i = 21). 

From all the pairs of adjacent eigenstates (e^e^+i), 
which give rise to the lowest frequency modes (coj i+i = 
(e^+i — 6i)/h) in the dynamics, we identify in Fig.jsjjthose 
which have the highest (in absolute value) amplitude 
^fi+i- For small gate voltages these are the eigenstates 
around the Fermi level but as V g increases the modes 
with the largest amplitudes tend to arise where avoided 
crossings of eigenstates occur. These are also the modes 
with frequencies close to the lowest possible frequency for 
that particular gate voltage. 

The actual highest amplitude at the last site, corre- 
sponding to the highlighted modes in Fig. [9j are de- 
picted as a function of V g in Fig. [lQ^ a), together with 

the quantity = ^J2i( A i^i-\-i) 2 • This represents 

the pessimistic estimate (considering all modes orthogo- 
nal in phase) of the total low-frequency amplitude corre- 
sponding to all pairs of adjacent eigenvalues. Such curves 
demonstrate the resonant nature of the dynamical site- 
occupation as a function of the gate voltage. Fig. [lQ^a) 
also shows that the low frequency charge excitation peaks 
for intermediate gate voltages (for this case V g ~ 2.57) 
can be related to the observed enhanced low-frequency 
coupling between the itinerant and localized spins for 
the intermediate V g range. Interestingly, the two A(V g ) 
curves coincide for very low and very high voltages, show- 
ing that in this case the entire dynamics is practically due 
to just one mode (marked in black in Fig. [9]). In the 
intermediate regime it is clear that there are more than 
one low-frequency modes which have a significant ampli- 
tude, although the one marked in Fig.[9]gives the leading 
contribution. 
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For a particular pair of adjacent eigenstates there are a 
few instances where the adiabatic eigenvalues come closer 
together as a function of the gate voltage. In Fig. 10 'b,c) 
we analyze the dependence of the frequency and the am- 
plitude corresponding to the pair of eigenvalues just be- 
low and above the Fermi level in the ground state (i = 15 
for the 30-atom chain at half- filling) . We find that the 
dynamical charge amplitude on the last site peaks (as 
absolute value) every time the frequency passes through 
a minimum. 



-l.l 



FIG. 9: (Color online) Adiabatic variation with the gate volt- 
age of a few eigenvalues around the Fermi level (we have used 
So = and | = 1 eV, hence Ep = at V g = 0). Marked with 
black lines are the pairs of eigenstates (i, 1) which give rise 
to the highest (in absolute value) dynamical amplitude A^ i+1 
at the last site N. 
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FIG. 10: (Color online) (a) The maximum single-mode am- 
plitude Aj| ax and an estimate for the total low-frequency am- 
plitude (see text for details) as a function of the applied 
gate voltage, (b) The frequency of the mode related to the 
pair of eigenstates (i = 15) around the Fermi level, adiabati- 
cally evolved with the gate voltag e and (c) its corresponding 
amplitude as defined in Eq. (A.l). 
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